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(Lecture Note for the SERC School held at Mehta Research Institute, India) Abstract. Low- 
dimensional quantum spin systems are interacting many body systems for which sev- 
eral rigorous results are known. Powerful techniques like the Bethe Ansatz provide ex- 
act knowledge of the ground state energy and low-lying excitation spectrum. A large 
number of compounds exists which can effectively be described as low-dimensional 
spin systems. Thus many of the rigorous results are also of experimental relevance. 
In this Lecture Note, an introduction is given to quantum spin systems and the rig- 
orous results known for such systems. Some general theorems and exactly-solvable 
models are discussed. An elementary introduction to the Bethe Ansatz technique is 
given. 

1. Introduction 

In Condensed Matter Physics, we study many body systems of different kinds. In a many body 
system, a large number of entities interact with each other and usually quantum effects are im- 
portant. The entities may be fermions (electrons in a metal, He'^ atoms ), bosons (He^'atoms) or 
localised spins on some lattice. In my set of lectures, we will discuss interacting spin systems in 
low dimensions. The spins are located at the different sites of a lattice. The most well-known 
model of interacting spins is the Heisenberg model with the Hamiltonian 

i7=^4.s:.^ (1) 

Si is the spin operator located at the lattice site i . Jij denotes the strength of the exchange 

interaction. The spin | Si |can have a magnitude 1/2, 1, 3/2, 2, etc. The lattice, at the sites 

of which the spins are located , is d-dimensional. Examples are a linear chain (d = 1), the square 
lattice (d = 2 ) and the cubic lattice (d = 3). Ladders have structures interpolating between the 
chain and the square lattice. A n-chain ladder consists of n chains coupled by rungs. The spin 
dimensionality n denotes the number of components of a spin. The Ising, XY and Heisenberg 
models correspond to n =1, 2 and 3 respectively. The generalised Ising and XY Hamiltonians are: 

Hising — JijSiSj (2) 

H,^=J2J^AS^S]+SfSy] (3) 

(ij) 

For the anisotropic XY model, J^^ ^ Jfj . Usually, the sites i and j are nearest-neighbours (n.n.s) on 
the lattice and the Jy 's have the same magnitude J for all the n.n. interactions. The Hamiltonian 
in (1) becomes 

H = JJ2S',.S', (4) 

There are, however, other possibilities. We mention some examples : 
(a) Alternating chain system : 

H = Y.[i-{-iys]s;.-^,+^ (5) 
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where < S <1. 

(b) Further-neighbour interactions : 

The well-known Majumdar-Ghosh chain is described by the Hamiltonian 

N N 

Hmg — J 

i=l i=l 

The Haldane-Shastry model[H has a Hamiltonian of the form 

(c) Jij 's can be random in sign and magnitude |^. 

(d) Quantum antiferromagnetic chain with incommensuration 

(e) Mixed-spin system: 

A good example is a chain of alternating spins 1 and 1/2 . 

(f) Fully anisotropic Heisenberg Hamiltonian : 
The Hamiltonian in Id and for n.n. interaction 



N 

HxYz = ^ [JxSfS^^i + JySfSf_^_^ + JzS^S^^i] (8) 

i=l 

The Ising and XY Hamiltonians are special cases of Hxyz ■ 

Once we know the Hamiltonian of a spin system, we want to determine the eigenvalues and 
the eigenfunctions of the system. This is, however, a difficult task as the number of eigenstates 
is huge . For a spin-1/2 system, the number is 2^ where N is the number of spins. We usually 
determine the ground state and the low-lying excited states in an approximate manner. Exact 
knowledge can be obtained only in a few cases. The low- lying excitation spectrum specifies the 
low-temperature thermodynamics and the response to weak external fields. The spectrum may 
have some universal features. In the case of electronic systems, there are two well-known universal 
classes: the Fermi liquids and the Luttinger liquids. We will learn about the universal features of 
spin excitation spectra as we go on. The spin excitation spectrum can be gapless or can have a gap. 
In the case of a gapless excitation spectrum, there is at least one momentum wave vector for which 
the excitation energy is zero. For a spectrum with gap, the lowest excitation is separated from 
the ground state by an energy gap A. The temperature dependence of thermodynamic quantities 
is determined by the nature of the excitation spectrum (with or without gap). Again, one can 
give examples from an electronic system. For a fermi liquid, the electronic excitation spectrum is 
gapless and the electronic specific heat C'n has a linear temperature dependence. A conventional 
superconductor is characterised by an energy gap A in the electronic excitation spectrum and the 
electronic specific heat c„ ~ exp—-^^ in the superconducting state. 

The partition function Z of a thermodynamic system is 

Z^^e"'^ (9) 

i 

where Ei 's are the energy levels and the sum is over all the states of the system. The free energy 

F = -ksTlnZ (10) 

The thermodynamic quantities of a magnetic system are: 
Magnetic energy : 

87 

Entropy Sh at constant field H : 
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Magnetization Mt at constant temperature 



OF 



MT = -{gjj] (13) 



Specific heat at constant field : 



Isothermal susceptibility 



fdM\ _ _ fd^\ 



XT=^ =-Ut72 (15) 



\dHjj, \dH^ J 



T 



Spontaneous magnetization Ms (H = 0) is known as the order parameter and has a non-zero value 
in the ordered phase of a magnetic system. Susceptibility, known as a response function of the 
system, is a measure of how the magnetization changes as the external field H is changed. The 
zero-field susceptibility is 

fdM 



At low temperatures, only the lowest energy levels are excited and a knowledge of the low-lying 
excitation spectrum enables one to determine the low temperature thermodynamics. 

Exchange interaction can give rise to magnetic order below a critical temperature. However, 
for some spin systems, there is no magnetic order even at T = 0. The three major types of 
magnetic order are ferromagnetism (FM), antiferromagnetism (AFM) and ferrimagnetism (Fi) . 
For FM order, the n.n. spins have a tendency towards parallel alignment, for AFM order, the n.n. 
spins tend to be antiparallel whereas for Fi order, the n.n. spins favour antiparallel alignment but 
the associated magnetic moments have different magnitudes. In all the three cases, the magnetic 
order survives upto a critical temperature Tc . The order parameter in the case of an AFM is 
the spontaneous sublattice magnetization. Another measure of the spontaneous order is in the 
two-spin correlation function. Long range order (LRO) exists in the magnetic system if 

liruR^^ (^(0).^(:^)\ ^ (17) 



where Tt denotes the spatial location of the spin. At T = 0, the expectation value is in the ground 
state and at T ^ 0, the expectation value is the usual thermodynamic average. 

The dynamical properties of a magnetic system are governed by the time-dependent pair cor- 
relation functions or their Fourier transforms. An important time-dependent correlation function 



G{R,t) = ^SR{t).So{0)) (18) 

Dynamical correlation function is the quantity measured in inelastic neutron scattering experi- 
ments. The differential scattering cross-section in such an experiment is given by 



dfldu) 



(xS'"'{-f,uj) (19) 



1 — > ~e! C^°° 

-5]e^9.i?/ rfte-*(5^(05o^(0)) (20) 

R J — CO 
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where ~^ and oj are the momentum wave vector and energy of the spin excitation. The spin 
components, (/i = x, y, z ), are perpendicular to . For a particular "g* , the peak in 
S^^{'q',uj) occurs at a value of a; which gives the excitation energy. At T = 0, 

S^^^'{t,Lo)=J2K^i^ + Eg-Ex) (21) 

A 

Eg {Ex ) is the energy of the ground (excited) state and 

= 2n\{G\ 5^( g>) I A) p (22) 

Quantities of experimental interest include the dynamical correlation function, various relax- 
ation functions and associated lineshapes. 

2. Ground and excited states of the FM and AFM Heisenberg 
Hamiltonian 

Consider the isotropic Heisenberg Hamiltonian 

H = J^'^i.S'j (23) 

{ij) 

where (ij) denotes a n.n. pair of spins. The sign of the exchange interaction determines the 
favourable alignment of the n.n. spins. .J > ( .J < 0) corresponds to AFM (FM) exchange 
interaction. To see how exchange interaction leads to magnetic order, treat the spins as classical 
vectors. Each n.n. spin pair has an interaction energy JS^ cos 9 where 9 is the angle between 
n.n. spin orientations. When J is < 0, the lowest energy is achieved when 9 = 0, i.e., the 
interacting spins are parallel. The ground state has all the spins parallel and the ground state 
energy Eg = —^^S'^ where z is the coordination number of the lattice. When J is > 0, the lowest 
energy is achieved when 9 = ir, i.e., the n.n. spins arc antiparallel. The ground state is the Neel 
state in which n.n. spins are antiparallel to each other. The ground state energy Eg = —^i^^S^. 

Magnetism, however, is a purely quantum phenomemon and the Hamiltonian (23) is to be 
treated quantum mechanically rather than classically. For simplicity, consider the case of a chain 
of spins of magnitude 1/2. Periodic boundary condition is assumed , i.e.,^jv+i = ■ 

N N 

H = J ^2 Si-'^i+i = J ^ 

i=l i=l 

where 



(24) 



Sf = S^±iSl (25) 

are the raising and lowering operators. For spins of magnitude ^ , has two possible values 1/2 
and —1/2 {h = 1). The corresponding states are denoted as \a) (up-spin) and |/3) (down-spin) 
respectively. The spin algebra is given by 





= \W) 


S^\I3) 


= 4"^^ 


S+\a) 


= 


5+1/3) 


= l«) 


S-\a) 


= m 




= 



(26) 
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It is easy to check that in the case of a FM, the classical ground state is still the quantum 
mechanical ground state with the same ground state energy. However, the classical AFM ground 
state (the Neel state) is not the quantum mechanical ground state. The determination of the exact 
AFM ground state is a tough many body problem and the solution can be obtained with the help 
of the Bethe Ansatz technique (Section 4). 

For a FM, the low-lying excitation spectrum can be determined exactly. Excitations are created 
by deviating spins from their ground state arrangement. Deviate a single spin in the FM ground 

state jaaao! aaad) . The number of such states is N as the deviated spin can be at any one of 

the N locations (sites). The actual eigenfunction is a linear combination of these states: 

JV 

= ^ a(m)V'(m) (27) 

m— 1 

where m denotes the location of the deviated spin and runs from 1 to N. The wave function i\)(m) 
has a down spin /? at the mth site. 

'0(m) = \aaaa /3aa....aa) (28) 

To solve the eigenvalue problem 

iJ* = (29) 
one has to determine the unknown coefficients a(m). 

H ^ a{m)tp{m) = -B X] «MV'M (30) 

m m 

Multiply by •^^ {I) on both sides to get 

Y,<m)m)\HW{m)) = EY,<rn)m)\^P{m)) 

m m 

m 

= Ea{l) (31) 

On the l.h.s., only those terms survive for which H \%l){m)) = \ip{l)) ■ Consider the FM version 
of the Heisenberg Hamiltonian in (24), i.e., replace J by -J with J> 0. The Hamiltonian H is a sum 
of two parts: 

H = H, + HxY (32) 

We consider the effect of each Hamiltonian on a n.n. spin pair. The following rules hold true. 
acting on a parallel (antiparallel) spin pair gives - J/4 (J/4) times the same spin pair. 

- J^f \aa) = \aa) , -JS^S^i = ^ \a(3) (33) 

HxY acting on a parallel spin pair gives zero and acting on an antiparallel spin pair interchanges 
the spins with the coefficient —J/2 

-^{StS-,, + S-St,,)\aa) = 
-^{StS-+, + SrS+,,)\a/3) = -^|/3a) 

~iStSr+i + S^Sl,)\l3a) = ~\a(3) (34) 



5 



In the wave function i/'(m) , there are {N~2) parallel and 2 antiparallel pairs of spins. Applying 
the rules of spin algebra, we get from (31), 



--(TV -2) + - X 2 



Ea{l) = 

ea{l) =1 [2a{l) - a{l + 1) - a{l - 1)] (35) 



where, e = E + ^ , is the energy of the excited state measured w.r.t. the ground state energy 
Eg = . The solution for a{l) in (35) is 

a(0 = e'''^ (36) 

and we get 

e = J(l - cos(fc)) (37) 

From the PBC, 

a(0 = a{l + N) (38) 

which leads to 

e'''^ = 1 (39) 

or, k= 25:a,A = 0, 1,2, , iV - 1 . 

The excited state consists of a deviated down spin propagating along the chain. The excitation 
is called a spin wave or a magnon. For r down spins (r = 2,3,4,... etc.), we have r magnon 
excitations. The r magnons can scatter against each other giving rise to a continuum of scattering 
states or they can form bound states. The exact bound state spectrum can be derived using the 
BA (Section 4). The case of the AFM Heisenberg Hamiltonian will be considered in this Section. 

The BA technique is applicable to quantum many body systems only in Id. However, not 
all Id quantum many body problems can be solved by the BA. For such problems as well as for 
models in higher dimensions, approximate techniques are used to gain knowledge of the ground 
state and the excitation spectrum H, ||. 

3. Theorems and rigorous results for Heisenberg AFMs 

A. Lieb-Mattis theorem [Q: For general spin and for all dimensions and also for a bipartite 
lattice, the entire eigenvalue spectrum satisfies the inequality 

EoiS) < EoiS + 1) 

where Eq{S) is the minimum energy corresponding to total spin S. The weak inequality becomes 
a strict inequality for a FM exchange coupling between spins of the same sublattice. The theorem 
is valid for any range of exchange coupling and the proof does not require PBC. The ground 
state of the S=l/2 Heisenberg AFM with an even number N of spins is a singlet according to the 
Lieb-Mattis theorem. 

B. Marshall's sign rule |§, §]: 

The rule specifies the structure of the ground state of a n.n. 8 = 1/2 Heisenberg Hamiltonian 
defined on a bipartite lattice. The rule can be generalised to spin S, n.n.n. FM interaction but 
not to n.n.n. AFM interaction. A bipartite lattice is a lattice which can be divided into two 
sublattices A and B such that the n.n. spins of a spin belonging to the A sublattice are located 
in the B sublattice and vice versa. Examples of such lattices are the linear chain, the square and 
the cubic lattices. According to the sign rule, the ground state ip has the form 



= E^A^i^) (40) 
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where is an Ising basis state. The coefficient has the form 

C^ = {-r^a^ (41) 

with real and >0 and is the number of up-spins in the A sublattice. 
Proof: We write the AFM Heisenberg Hamiltonian in the form 



<ij> 

= 2J ^ {S!S] -\) + JY1 [S^Sr + SiS+] (42) 

<ij> <ij> 

= H,+HxY (43) 

Let rn^ be the z-component of the total spin and the number of antiparallel bonds in the spin 
state l/i). Hz acting on a parallel spin pair gives zero and acting on an antiparallel spin pair gives 
—J. HxY acting on a parallel spin pair gives zero and interchanges the spins in the case of an 
antiparallel spin pair. In the latter case, and are unchanged and ^ ± 1 . 

iJ|^) = J^C4-g^|M) + ^k)} (44) 
where \v) differs from by one interchange of antiparallel n.n.s. There are such states. 



P {^\H\^) 



I c, p 



Let 



1^ 1^ 



Also, C^C* = a^a*. From (45), 



(45) 



C, = (-ly-a, (46) 
where is arbitrary. Since p,^ = ± 1, from (46), 



(47) 



Now, since yit, v are dummy indices 

EE"X = ^EEM-+«>] (49) 

^ v pi/ 

If a^^a;/ are real and positive, —^\ |, ^| Ui, |and 

a^a^ =1 II fli/ I (50) 
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If a^a^ are complex, the term in (49) is 

\{\ I e*^- I a, \ e-"^-^+ \ a,, \ e"''"^- | a, \ e'^-} =| || a, | cos(0^ - 0,) 



< I Op 1 1 fli. I 

The energy in (45) can be reduced by choosing a^'s to be real and >0 (proved). 
C. Lieb, Schultz and Mattis theorem|§, ||]: 

A half-integer S spin chain described by a reasonably local Hamiltonian respecting translational 
and rotational symmetry either has gapless excitation spectrum or has degenerate ground states, 
corresponding to spontaneously broken translational symmetry. 

Proof: We consider the Heisenberg AFM Hamiltonian in (24). A rigorous proof of uniqueness 
(non-degeneracy) of the ground state |V'o) exists in this case. We wish to prove that there is a low 
energy excitation of 0(1/L), where L is the length of the chain. The proof consists of two steps : 
(i) Construct a state iV'i) which has low energy, i.e.. 



(^i|(i/-£;o)|Vo) = 0(1/L) (51) 

(ii) Show that is orthogonal to the ground state |V'o)- The state is constructed by applying 
an unitary transformation on which gives rise to a slowly varying rotation of the spins about 
the z-axis. 



U = exp 



(52) 



The operator U rotates the nth spin by the small angle about the z-axis. The spin coordinates 
are rotated by 2n about the z-axis between the first and last sites. Note the spin identities 



S-e 



Now, 



(53) 



{'iPi\{H - Ea)\^i) = {MU-\H - Eo)U (54) 
The Hamiltonian H in (24) consists of two parts and Hxy ■ Elz commutes with U. Therefore, 



= {M H iVo) - {MHxY iV'o) + {M U-^HxyU l^n) 

= Eo-{\{StS;+^)+h.c.} + {MU-^HxYU\i^o) (55) 

In the last term, making use of the spin identities, we get 



(56) 



So, from (54) - (56), 



(57) 
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Since the ground state is unique it must be isotropic implying 



(58) 



Thus, 

{i^i\iH-Eo)\^Pi) - 
Since S^S^_^i is a bounded operator, 



cos{ — ) - 1 



{^i\{H-Eo)\i^i) = 0{l/L) 



(59) 



(60) 



In the limit L ^ oo , the r.h.s. of (60)^ 0. 

Constructing a low-energy state, however, proves nothing. The state may become equal to the 
ground state as L ^ oo . To complete the proof, we have to show that this does not happen. This 
can be done by showing that are \ipo) are orthogonal to each other. We will now show that 
\ipi) has momentum tt relative to the ground state |V'o) which , being unique, is also a momentum 
eigenstate. Thus the two states are orthogonal to each other. Let T be the translation operator 
which produces a translation by one site. 



(61) 



TUT- 



L-l 



= exp 



= U exp 



.2tt 



(62) 



The ground state has total spin zero since it is unique and so the first exponential has the value 
1. Thus, 



T|^i) - TU\^o) ^TUT-'T\^Po) 

= e^^'^'^if/TlV'o) (63) 

The exponential is -1 since Sf has eigenvalues ±i . So, 

TlV'i) = -UT\i,o)=e'^^+''"'>U\^o} 

= e''("+^«) l^i) (64) 

where kg is the momentum wave vector of the ground state. The state has momentum 
wave vector tt relative to the ground state \4>o) and is thus orthogonal to it. This completes 
the proof of the LSM theorem for the S=l/2 Heisenberg AFM Hamiltonian in Id. The proof 
extends immediately to arbitrary half-odd integer S but fails in the integer S case. For the latter, 
the exponential in (63) is +1 and it cannot be proved that \ipi) and \ipo) are orthogonal to each 
other. If the translational symmetry is spontaneously broken in the ground state, the ground 
state develops degeneracy (Eq.(60) is still true but is another ground state as L ^ oo ). The 
excitation spectrum has a gap in this case. 

In the presence of a magnetic field, the S=l/2 HAFM chain has gapless excitation spectrum 
from zero field upto a saturation field, at which the ground state becomes the fully polarized 
ferromagnetic ground state. For integer S, Haldane made a conjecture that the excitation spectrum 
has a gap. This conjecture has been verified both theoretically and experimentally In the 
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presence of a magnetic field, the gap persists up to a critical field, equal to the gap. The spectrum 
is gapless from the critical field to a saturation field. The Lieb-Scultz-Mattis (LSM) theorem can 
be extended to the case of an applied magnetic field ||l^. In this case, a striking phenomenon 
analogous to the quantum Hall effect occurs. The phenomenon is that of magnetization plateaus. 
The content of the extended theorem is : translationally invariant spin chains in an applied field 
can have a gapped excitation spectrum, without breaking translational symmetry, only when the 
magnetization/site m obeys the relation 



S ~ m — integer (65) 

where S is the magnitude of the spin. The proof is an easy extension of that of the LSM theorem. 
The gapped phases correspond to magnetization plateaus in the m vs. H curve at the quantized 
values of m which satisfy (65). Whenever there is a gap in the spin excitation spectrum, it is 
obvious that the magnetization cannot change in changing external field. Fractional quantization 
can also occur, if accompanied by (explicit or spontaneous ) breaking of the translational symmetry. 
In this case, the plateau condition is given by 



n[S — to) = integer (66) 

where n is the period of the ground state. Hida jl^ has considered a S=l/2 HAFM chain with 
period 3 exchange coupling. A plateau in the magnetization curve occurs at m=l/6 (1/3 of full 
magnetization). In this case, n = 3, S = 1/2 and m = 1/6 and the quantization condition in (66) 
is obeyed. 

D. Mermin- Wagner's theorem : 

There cannot be any AFM LRO at finite T in dimensions d = 1 and 2. AFM LRO can, however, 
exist in the ground states of spin models in d = 2. LRO exists in the ground state of the 3d 
HAFM model for spin S> 1/2 |14|. At finite T, the LRO persists upto a critical temperature Tc. 
For square ||l5j and hexagonal |16| lattices, LRO exists in the ground state for S> 1. The above 
results are based on rigorous proofs. No such proof exists as yet for S=l/2, d =2 ( this case is of 
interest because the Cu02 plane of the high-Tc cuprate systems is a S = 1/2 2d AFM ). 
3. Possible ground states and excitation spectra 

The Neel state is not the ground state of the quantum Heisenberg Hamiltonian. the Neel 
state, we know, has perfect AFM LRO. It is still possible that the true ground state of a quantum 
Hamiltonian has Neel-type order, i.e., non-zero sublattice magnetization. The ground state is 
essentially the Neel state with quantum fiuctuations mixed in. This picture does indeed describe 
the ground state of the HAFM in d > 2. The ground state of the HAFM in Id has no LRO and is 
an example of a spin liquid . Even in higher dimensions, one can have disordered ground states due 
to what is known as frustration. Frustration occurs when the interaction energies associated with 
all the spin pairs cannot be simultaneously minimized. This may occur due to lattice topology 
or due to the presence of further-neighbour interactions. The triangular jl^, kagome[]l8| and 
pyrochlore lattices ||l^ are well-known examples of frustrated lattices. Another recent example is 
the pentagonal lattice |Q. In these lattices, the number of bonds in each elementary plaquette 
is odd. We illustrate frustration by considering Ising spins on a triangular plaquette. Each Ising 
spin, sitting at a vertex of the triangle, can be either up or down. AFM exchange interaction 
energy is minimized if the n.n. spin pairs are antiparallel. As Figure 1 shows, not all the three 
n.n. spin pairs can be simultaneously antiparallel. The parallel spin pair may be located on any 
one of the bonds of the triangle. In fact, the ground state of the Ising AFM on the triangular 
lattice is highly degenerate. The entropy/spin has a non-zero value in the ground state and there 
is no ordering of the spins. If there is a mixture of FM and AFM interactions, then a spin model 
defined on an unfrustrated lattice like the square lattice (each elementary plaquette, a square, has 
four bonds) can also exhibit frustration, if the number of n.n. AFM interactions in a plaquette is 
odd. Frustration also occurs if both n.n. as well as n.n.n. interactions are present. Frustration 
favours disordered ground states but does not rule out AFM LRO in the ground state. The Ising 
AFM on the triangular lattice has a disordered ground state but various calculations suggest that 
the HAFM Hamiltonian on the triangular lattice has AFM LRO [|l). 
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Various types of disordered ground states are possible: dimer, resonating valence bond (RVB), 
twisted, chiral and strip or collinear states. The lowest energy configuration of two spins interacting 
via the HAFM exchange interaction is a spin singlet with energy -3J/4. The singlet is often termed 
a valence bond (VB) or dimer (an object which occupies two sites). A pictorial representation 
of a VB between spins at sites i and j is a solid line joining the sites. The square lattice with 
N sites can be covered by N/2 VBs in various possible ways. Each possibility corresponds to a 
VB state. A RVB state is a coherent Hnear superposition of VB states (Figure 2). The 2d S 
= 1/2 AFM CaV40g has a lattice structure corresponding to a 1/5-depleted square lattice. The 
lattice consists of plaquettes of four spins connected by bonds. The ground state of the AFM is 
the Plaquette RVB (PRVB) state [22]. In this state, the spin configuration in each plaquette is a 
linear combination of two VB configurations. In one of these the two singlet bonds are horizontal 
and in the other they are vertical. Bose and Ghosh [23] have constructed a spin model for which 
the PRVB state is the exact ground state. The triplet excitation is separated from the ground 
state by an energy gap and singlet excitations fall within the gap. 

In a columnar dimer state, the dimers or VBs are arranged in columns. A description of the 
other disordered ground states can be obtained from Ref. [24]. The disordered ground states do 
not have conventional LRO but are not simple paramagnets. One can define suitable operators, 
the expectation values of which in the appropriate ground states are non-zero. The states, with 
different types of order, are described as quantum paramagnets. There is a huge literature on 
AFM spin models with disordered ground states. This is mainly because of the relevance of such 
studies for high temperature superconductivity. The common structural ingredient of the high-Tc 
cuprate systems is the Cu02 plane. It looks like a square lattice with Cm^"*" ions located at the 
lattice sites and oxygen ions sitting on the bonds in between. The Cv?^ ions carry spin-1/2 and 
the n.n. spins interact via the HAFM exchange interaction Hamiltonian. In the undoped state, 
the cuprates exhibit AFM LRO below the Neel temperature. On doping the systems with a few 
percent of holes, the AFM LRO is rapidly destroyed leaving behind a spin disordered state. It 
is important to understand the nature of the spin disordered state as this state serves as the 
reference state into which holes are introduced on doping. The holes are responsible for charge 
transport in the metallic state and forms bound pairs in the superconducting state. We now 
discuss some spin models for which the ground state properties (energy and wave function as well 
as correlation functions) can be calculated exactly. The first model is described by the Majumdar- 
Ghosh Hamiltonian [|l|, ^. This is defined in Id for spins of magnitude 1/2. The Hamiltonian is 
AFM in nature and includes both n.n. as well as n.n.n. interactions. The strength of the latter is 
half that of the former. 

N N 

i^A/G = J^^..^.+i+^/2^^..^,+2 (67) 

i=l 2=1 

The exact ground state of Hmg is doubly degenerate and the states are 



01 = [12] [34] [iV-liV] 

02 = [23] [45] [Nl] (68) 

where [Im] denotes a singlet spin configuration for spins located at the sites 1 and m. Also, PBC is 
assumed. One finds that translational symmetry is broken in the ground state. The proof that 
and 02 are the exact ground states can be obtained by the method of 'divide and conquer'. One 
can verify that 0i and 02 are exact eigenstates of Hmg by applying the spin identity 

S^„.(S^i+ S^„0[H -0 (69) 
Let El be the energy of 0i and 02 . Let Eg be the exact ground state energy. Then 

Eg < El (70) 
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One divides the Hamiltonian H into sub-Hamiltonians, Hi 's, such that H = J2i^i- ^^"^^ be 
exactly diagonalised and let -Ejo be the ground state energy. Let ■^g be the exact ground state 
wave function. By variational theory, 

i i 

From (70) and (71), one gets 

Y,Eio<Eg<Ei (72) 

i 

If one can show that E^ = Ei , then E\ is the exact ground state energy. For the MG-chain, 
the sub-Hamiltonian Hi is 

Hi = J/20i.'^ + ~§i+2.~§i) (73) 

There are N such sub-Hamiltonians. One can easily verify that £"20 = —SJ/8 and Ei = — 
(-3 J/4 is the energy of a singlet and there are N/2 VBs in 4>i and (^2 )• From (72), one finds that 
the lower and upper bounds of Eg are equal and hence and (p2 are the exact ground states with 
energy Ei = —3JN/8 . There is no LRO in the two-spin correlation function in the ground state: 

K\^,J) = {S!S^) = y,,-l^^i_J^,, (74) 
The four-spin correlation function has off-diagonal LRO. 



= K'^{ij)K'^{lm) + -^(5|i_j|,i5|i_„|,i 

X earp(i7r(^ - ^-)) (75) 
Let T be the translation operator for unit displacement. Then 

T(j)i = (l>2,T(j)2 = (j)! (76) 

The states 

</,+ = -^(01+</.2),</)- = -^(</.l-<A2) (77) 

correspond to momentum wave vectors k = and k = tt . The excitation spectrum is not exactly 
known. Shastry and Sutherland [25] have derived the excitation spectrum in the basis of 'defect' 
states. A defect state has the wave function 



V'(p,m) = ....[2p-3,2p-2]a2p-i[2p,2p+l] 

....[2m -2,2m- l]a2m[2m + l,2m + 2] (78) 

where the defects {a2p-i and a2m ) separate two ground states. The two defects are up-spins and 
the total spin of the state is 1. Similarly, the defect spins can be in a singlet spin configuration so 
that the total spin of the state is zero. Because of PBC, the defects occur in pairs. A variational 
state can be constructed by taking a Hnear combination of the defect states. The excitation 
spectrum consists of a continuum with a lower edge at J {5/2 — 2 | cosk |). A bound state of 
the two defects can occur in a restricted region of momentum wave vectors. The MG chain has 
been studied for general values a J of the n.n.n. interaction [26]. The ground state is known 
exactly only at the MG point a = 1/2 .The excitation spectrum is gapless for < a < Ucr (— 
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0.2411). Generalizations of the MG model to two dimensions exist [27,28]. The Shastry-Sutherland 
model [27] is defined on a square lattice and includes diagonal interactions as shown in Figure 3. 
The n.n. and diagonal exchange interactions are of strength Ji and J2 respectively. For ^ below a 
critical value ^ 0.7, the exact ground state consists of singlets along the diagonals. At the critical 
point, the ground state changes from the gapful disordered state to the AFM ordered gapless 
state. The compound SrCu2{B03)2 is well-described by the Shastry-Sutherland model [29]. Bose 
and Mitra [28] have constructed a Ji — J2 — J3 — J4 — ^5 spin-1/2 model on the square lattice. 
Ji,J2,J3,J4 and J5 are the strengths of the n.n., diagonal , n.n.n., knight 's-move-distance-away 
and further-neighbour -diagonal exchange interactions (Figure 4) . The four columnar dimer states 
(Figure 5) have been found to be the exact eigenstates of the spin Hamiltonian for the ratio of 
interaction strengths 

Ji : J2 : J3 : J4 : ^5 = 1 : 1 : 1/2 : 1/2 : 1/4 (79) 

It has not been possible as yet to prove that the four columnar dimer states are also the ground 
states. Using the method of 'divide and conquer', one can only prove that a single dimer state is 
the exact ground state with the dimer bonds of strength 7J. The strengths of the other exchange 
interactions are as specified in (79). 

In Section 2, we have discussed the LSM theorem, the proof of which fails for integer spin chains. 
Haldane [30,10] in 1983 made the conjecture, based on a mapping of the HAFM Hamiltonian, in 
the long wavelength limit, onto the nonlinear a model, that integer-spin HAFM chains have a 
gap in the excitation spectrum. The conjecture has now been verified both theoretically and 
experimentally [31]. In 1987, Affleck, Kennedy, Lieb and Tasaki (AKLT) constructed a spin-1 
model in Id for which the ground state could be determined exactly [32]. Consider a Id lattice, 
each site of which is occupied by a spin-1. Each such spin can be considered to be a symmetric 
combination of two spin-l/2's. Thus, one can write down 



= \++),S^ = +l 
V- = \—),S^ = -l 

V+- = -^(|+-) + |-+)),5^ = 

= V'+- (80) 

where ('— ') denotes an up (down) spin. 

AKLT constructed a valence bond soHd (VBS) state in the following manner. In this state, each 
spin-1/2 component of a spin-1 forms a singlet (valence bond) with a spin-1/2 at a neighbouring 
site. Let e"'^ = +ot — ) be the antisymmetric tensor : 

e++ = e~ =0,e+- = -e-+ = 1 (81) 

A singlet spin configuration can be expressed as -^s"'' I"/?) , summation over repeated indices 
being implied. The VBS wave function (with PBC) can be written as 



i^a,0,e'^'"'+' ^a^p^e^""' (82) 

I'fpVBs) is a linear superposition of all configurations in which each 5^ = -M is followed by a 
= —1 with an arbitrary number of = spins in between and vice versa. If one leaves out 
the zero's, one gets a Neel-type of order. One can define a non-local string operator 

ag. = -Srexpiin ^ 5^)5^, (a = x, y, z) (83) 
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and the order parameter 



0?trms = ^«m|,_j|^oo (ffg-) (84) 

The VBS state has no conventional LRO but is characterised by a non-zero value 4/9 of O^trmg- 
After constructing the VBS state, AKLT determined the Hamiltonian for which the VBS state is 
the exact ground state. The Hamiltonian is 

Haklt = P2(Sl + ^i+i) (85) 

i 

where P2 is the projection operator onto spin 2 for a pair of n.n. spins. The presence of a VB 
between each neighbouring pair imphes that the total spin of each pair cannot be 2 (after two of 
the S=l/2 variables form a singlet, the remaining S=l/2's could form either a triplet or a singlet). 
Thus, Haklt acting on \^pvBs) gives zero. Since Haklt is a sum over projection operators, the 
lowest possible eigenvalue is zero. Hence, IVvss) is the ground state of Haklt with eigenvalue 
zero. The AKLT ground state (the VBS state) is spin-disordered and the two-spin correlation 
function has an exponential decay.The total spin of two spin-l's is 2, 1 or 0. The projection 
operator onto spin j for a pair of n.n. spins has the general form 



[l{l 



[1(1 + 1) + 

where = 5^ + S^j+i. For the AKLT model, j =2 and 1 = 1,0. From (85) and (86), 

(87) 



Haklt = ^ 



The method of construction of the AKLT Hamiltonian can be extended to higher spins and to 
dimensions d >1. The MG Hamiltonian (apart from a numerical prefactor and a constant term) 
can be written as 



H = J2P3/20^ + ^^+l+^^+2) 



The S=l HAFM and the AKLT chains are in the same Haldane phase, characterised by a gap in 
the excitation spectrum. The physical picture provided by the VBS ground state of the AKLT 
Hamiltonian holds true for real systems [33]. The excitation spectrum of Haklt cannot be 
determined exactly. Arovas et al [34] have proposed a trial wave function 



N 



\k) = N-'^Y.^'''^J IV'VBs) ,M = (89) 

and obtained 



^ {k\HvBs\k) ^ 25 + irnMk) 
^ ' <k\k> 27 ^ ' 

The gap in the excitation spectrum A = ^ at k = tt. Another equivalent way of creating 
excitations is to replace a singlet bond by a triplet spin configuration [35]. 

We end this Section with a brief mention of spin ladders. The subject of spin ladders is a 
rapidly growing area of research [36,37,38]. Ref. [38] gives an exhaustive overview of real systems 
with ladder-like structure and the results of various experiments on ladder compounds. A spin 
ladder consists of n chains ( n = 2,3,4, ...etc.) coupled by rungs. The simplest spin ladder (Figure 
6) corresponds to n = 2. Each site of the ladder is occupied by a spin- 1/2. Let .J and Jr be 
the strengths of the intrachain n.n. and rung exchange interactions. When Jr = , the ladder 
decouples into two HAFM chains for which the excitation spectrum is known to be gapless. A gap 
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opens up in the excitation spectrum even for an infinitesimal value of Jr and persists for all values 
of Jr . When Jr is >> J, the ground state predominantly consists of singlets along the rungs 
of the ladder. The ground state is spin-disordered. A triplet excitation is created by replacing a 
rung singlet by a triplet and letting it propagate. The n-chain spin ladder exhibits the 'odd-even' 
effect. The excitation spectrum is gapless (with gap) when n is odd (even). When n is odd, the 
spin-spin correlation function has a power-law decay as in the case of the HAFM S = 1/2 chain. 
For n even, the ground state correlation function has an exponential decay. For n sufBciently large, 
the square lattice limit is reached with AFM LRO in the ground state and a gapless excitation 
spectrum. Bose and Gayen [39] have constructed a two-chain spin ladder model which includes 
diagonal exchange interactions of strength J. When Jr >> 2 J, the exact ground state consists of 
singlets along the rungs. Later, Xian [39] has shown that the same exact ground state is obtained 
for ^ > ~ 1.40148. Ghosh and Bose [40] have generahsed the two-chain ladder model to an 
n-chain model for which the 'odd-even' effect can be exactly demonstrated. 

The spin ladder can be doped with holes. Two holes can form a bound state as in the case of 
high-Tc cuprate systems. The binding of two holes in a t-J ladder model has been shown through 
exact, analytic calculations [41]. A ladder compound has been discovered which shows supercon- 
ductivity (maximum Tc ~ 12K ) on doping with holes and under pressure. In the superconducting 
state, the holes form bound pairs. Again, the analogy with the cuprate systems is strong. For 
detailed information on undoped and doped spin ladders, see the Refs. [36,38]. 
4. The Bethe Ansatz 

The Bethe Ansatz (BA) was formulated by Bethe in 1931 [42] and describes a wave function 
with a particular kind of structure. Bethe considered a well-known model in magnetism, the spin- 
1/2 Heisenberg linear chain in which only n.n. spins interact. The wave function of the interacting 
many body system has the BA form. For ferromagnetic (FM) interactions, Bethe starting with 
the BA wave function derived the energy dispersion relations of spin wave (magnon) bound states 
exactly. Hulthen [43] used the BA method to derive an exact expression for the ground state 
energy of the antiferromagnetic (AFM) spin-1/2 Hamiltonian. In later years, it was realized that 
the method has a wider applicability and is not confined to magnetic spin chains. In fact, many 
of the exact solutions of interacting many body systems are but variations or generahzations of 
Bethe's method. Examples include the Fermi and Bose gas models in which particles on a line 
interact through delta function potentials [44], the Hubbard model in Id [45], Id plasma which 
crystallizes as a Wigner solid [46], the Lai-Sutherland model [4 7] which includes the Hubbard model 
and a dilute magnetic model as special cases, the Kondo model in Id [48], the single impurity 
Anderson model in Id [49], the supersymmetric t-J model (J =2t) [50] etc. In the case of quantum 
models, the BA method is applicable only to Id models. The BA method has also been applied 
to derive exact results for classical lattice statistical models in 2d. In the following we show how 
the BA method works by considering the spin-1/2 linear chain Heisenberg model mentioned in the 
beginning. The Hamiltonian for the linear chain is given by 



describing the interaction of spins located at n.n. lattice sites. The total number of spins is given 
by N. The magnitude of the spin at each site is given by 1/2. Periodic boundary condition (PBC) 
is assumed so that N + 1 =1. We consider the interaction to be ferromagnetic so that the exchange 
constant J is > 0. The Hamiltonian describes magnetic insulators in which magnetic moments are 
localized on well-separated atoms. The Hamiltonian (91) can be written as 



N 




(91) 



1=1 




(92) 



There are two conserved quantities: 
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N 

= Y,S!,[S\H]=0 

i=l 

N 

~§ = Y.^i,[^\H] = (93) 

i=l 

Since 5^ is a conserved quantity, there is no mixing of subspaces corresponding to the different 
values of , thus making it easier to solve the eigenvalue problem. Define the spin deviation 
operator 

n=^ + S^ (94) 
The FM ground state V'g corresponds to n =0, i.e., = with all the spins pointing downwards: 

= mm (95) 

The ground state energy is Eg = — J-j . 

For n = r, we have r up -spins at the sites m\,m2, rrir. The eigenfunction ■i/' of H in the n 

= r subspace is a linear combination of the ^Cr functions ■^(mi,m2, jfrir) : 

ij) = ^a(mi,m2, ,mr)tp{mi,m2, ,mr),{Tn} = {mi,m2, ,mr) (96) 

{m} 

The summation is over all the mj 's running over 1 to N subject to the condition 



mi < m2 < < nir (97) 

to avoid overcounting of states. We now consider the operation of the Hamiltonian H given in (92) 
on V'- The rules of the spin algebra are: 

(a) if two neighbouring spins are parallel, the term S^S~ + S~ in H gives no contribution. 
The S^S^ term on the other hand gives — -j times the same state ■0({m}) , 

(b) for an antiparallel n.n. spin pair, the term + interchanges the spins in ?/>({m}) 
with a multiplicative constant — ^ . The S^S^ term gives a contribution j multiphed by the same 
state tl>{{m}) . The eigenvalue equation is given by 

Hi; = E^p = E^a{{m})^p{{m}) (98) 

{m} 

Take the scalar product with a particular il>*{mi,m2, ,TOr) and use orthogonality properties 

to get 

Ea{mi,m2, , m^) = -| ^ a{{m'}) + ^N'a{{m}) - '-^{N - N')a{{m}) (99) 

{m'} 

N' is the total number of antiparallel spin pairs so that {N - TV' ) is the number of parallel spin 

pairs. The sum is over all N' distributions {m[,m'2, ,m'j.) which arise from a n.n. exchange of 

antiparallel spins in (mi,m2, ,'fnr). Eq. (99) can further be written as 



[E + JN)a{mi,m2, ,"ir) = —4 [a("T-'i,W2, ,m^) — a(TOi,TO2, ,TOr)] (100) 



{m'} 



Or 
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2eo(mi,m2, ,mr) = — ^ [a(TO'i,m2, ,m'^) — a{mi,m2, ,mr)] (101) 

where 

e=j{E+^J) (102) 

is the energy, in units of J, of the excited states w.r.t. that of the ground state energy Eg. Next, 
one considers the PBC from which one gets 

a(mi, ....,mi, ,mr) = a(mi, ,mr,mi + N) (103) 

The spin waves of a FM system are obtained for n =1. 

^(m) = ^a{m)il){m) (104) 

m 

2ea{m) = - [a{m + 1) + a(m - 1) - 2a(m)] (105) 
The solution of (105) is given by 

a(m) = e''=™ 

and 

e = 1 - cosk (106) 

From PBC, the allowed values of k are 

fc=^,A = 0,l,...,iV-l (107) 

There are N states. Consider now the case n =2, i.e., there are two spin deviations. We have to 
take into account two cases. If the two spin deviations are not neighbours, one gets 

2ea(mi,m2) = — a(mi + l.m2) — a(mi — l,m2) — a(mi,m2 + 1) 

— a(mi, m2 — 1) + 4a(mi, m2) (108) 

If the spin deviations are neighbours, one gets 

2ea{mi, mi + 1) = — a(mi — 1, mi + 1) — o(mi, mi + 2) + 2a(mi, mi + 1) (109) 
Eq. (108) is satisfied by the Ansatz 

a(mi, ma) = Cie^C^i^i+'^^ms) ^ (J^^^{k2m^+k^m2) f^^^Q^ 

Then 

e = 1 — coski + 1 — cosk2 (HI) 

Ci, C2, fci and k2 are to be determined. Eq. (109) can be satisfied by Eq. (110), if the coefficients 
Ci and C2 are chosen to make 

a(mi,mi) + a(mi + l,mi + 1) — 2a(mi,mi + 1) = (112) 

For spin-1/2 particles, the amplitudes a(m,m) have no physical significance but are defined by Eq. 
(112). Put Eq. (110) into Eq.(112) to get 
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Ci sin^{ki ~ k2) + i{cos^{ki + — cos^{ki — k2)} (113) 
C2 sin^{ki — k2) — i{cos^{ki + ^2) — cos^{ki — k2)} 

Let Ci = e'i and C2 = e"*! Then 

2cot- = cot— - cot— (114) 
2 2 2 ^ ' 

Also, PBC gives 

a{mi,m2) = a{m2,mi + N) (115) 
This is in accordance with the ordering of mi 's given in Eq. (97). From the PBC, one gets 

Nki -(j) = 27rAi 
Nk2 + (f) = 2n\2 

Ai,A2 = 0,1,2, ,{N-1) (116) 

We have three equations (114) and (116) for the three quantities (f)\,ki,k2 ■ The sum of k\ and 
k2 is a constant of motion by translational symmetry. 

277 

fc = fci + fc2 = -^(Ai+A2) (117) 

Since Ai and A2 can be interchanged without affecting the solution, we choose Ai < A2. For a 
given Ai and A2, the solutions are completely determined by (114) and (116). 

The following three cases are to be considered for determining the dependence of the wave 
numbers fci and ^2 on Ai and A2 : 
Case I: A2 > Ai + 2 
Case II: Ai = A2 = -J , A even 
Case III: Ai = A2 - 1 = A odd 

In Case I, for a given value of A2, Ai can take on the values Ai = 0,1,2,....,A2 -2. The total number 
of solutions is given by 

7V-1 

^ (A2 - 1) =^-^ C2 (118) 

A2=2 

In this case, all the wave numbers are real and for specified values of Ai and A2 , fci, ^2 and (f) are 
determined uniquely from (114) and (116). The energy eigenvalue e/ is given by (111) with A;i,A;2 
real. In Case II, from (116), 



Nki -(j) = wX 

Nk2 + 4> = ttA (119) 



We look for solutions of the form 



fci = u + iv 
k2 = u — iv 

(j) = ii + ix (120) 

From (119) and (120), we find that tjj = and x = . If u is non-zero, then for N ^ 00, x also 
— > 00. So 
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In a first approximation, 



cot- 



sinip — isinhx 
coshx — cosip 
-z(l + 2e-x) 



Hence, 



2cot- 
2 

-2i 
or, sinhv 

— V 

or, e 



M k2 

cot cot 

2 2 

sinu — isinhv sinu + isinhv 



coshv — cosu 
coshv — cosu 

cosu 



coshv — cosu 



(121) 



2 — cos{u + iv) — cos{u — iv) 

2 — cosu{cosu H ^— ) 

cosu 

^(1 — cos2u) — ~ cosk) 



(122) 



From (121), cosu>0, so — 7r/2 < m < 7r/2 and k can be chosen to be in tiie range — tt < fc < tt. 
Let us specify the energy obtained in (122), for complex values of fci and k2 , as e/j . Let us now 
compare e/ and ejj . From (111), 



e/ = 1 — coski + 1 — cosk2 

= 2 — 2cos^{ki — k2)cos^{ki + k2) 
1 k 

= 2- 2cos-{ki - k2)cos- (123) 

The value of k lies in the range — tt < k < tt . For each value of k, ki — k2 lies in the range 
< (A:i — k2) < 2tt . Thus if one plots e/ against k, one gets a continuum of scattering states 
(Figure 7). 

Now, the minimum value of e/ is 

e™™ = 2- 2cos^ (124) 



From (122), 



Thus, 



e/7 = l-cos2^ (125) 



i(l + co4)<l (126) 



e 



mm 



2' 2' 



We find that e//<e™*" for k 7^0. For k= 0, the bound state is degenerate with an eigenstate 
belonging to the continuum of scattering states and having the energy eigenvalue zero. Thus there 
are N-1 bound state solutions. 

We now consider the wave functions in Cases I and IL Fos Case I, ki, k2 and (j) are real and 

a{mi,m2) = e'(*^i'"i+'=2™2+<)!>/2) ^ ^i{k2rm+k^m2-4>/2) ^^37) 
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The state describes the scattering of spin waves. The 'phase shift' cf) results from the mutual 
interaction between the spin waves. Since spin operators at two different sites commute, the 
excitations behave like bosons and the amplitude a(jni,m2) is properly symmetrized to reflect the 
bosonic nature of the spin waves. In Case II, from (120) and with ip = and x=Nv, we get 



a(mi,m2) = g'''t'/'^g''{kimi+k2m2) _^ ^-i(f>/2^i{k2mi+kim2) 
_ gife/2(TOi+r 



^—Nv/2—v{mi—m2) _|_ ^Nv/2+v{mi—m2) 



y - (m2 - mi) 



(128) 



If we normahse these states, we see that the two reversed spins tend to be localised at n.n. positions; 
I a(mi, m2) |is a maximum for TO2 = mi + 1 and has an exponential decay for m2 > mi + 1 . The 
width of the bound state is given by v and u is the velocity of the centre of mass of the bound 
complex in the chain. The Case II solutions represent the bound states of two spin waves. 

Finally, we consider Case III, i.e., Ai = (A — l)/2 and A2 = (A + l)/2 (A odd). From (116), we 

get 



Nki-(f) = 7r(A-l) 
Nk2 + (l> = 7r(A + l) 

A solution of (129) is fci = ^2 and ^ = tt. Then from (127), 



a(mi, m2) 



ife/2(mi+m2) 



[exp{i'!T/2) + exp{—m/2] = 



(129) 



(130) 



For complex fci,A;2 and (j), tp = n and x = Nv from (120). For N ^ 00 , we obtain the same 
solutions as we find in Case II. 

For the general case n = r, Bethe proposed a form for the amplitude a({m}) (the celebrated 
BA) which is a generalisation of the amplitude in the n =2 case. 



a{mi,m2, ,Tnr) = exp 

p 



3=1 



I > kpjruj 



1 ■ 



PI 



3<l 



(131) 



The sum over P denotes a sum over all the r! permutations of the integers l,2,...,r, Pj is the 
image of j under the permutation P. Each term in (131) has r plane waves scattering against one 
another. For each pair of plane waves, there is a phase shift <f>j^i . The sum over permutations is 
in accordance with the bosonic nature of the waves, the spin waves, propagating along the chain. 
For r spin deviations several possibilities occur: none of the spin deviations occur at n.n. sites, 
only two of the spin deviations occur at n.n. sites, more than a single pair of spin deviations occur 
at n.n. sites, three spin deviations occur at n.n. sites and so on. However, only the first two 
possibilities need to be considered to solve the eigenvalue problem, the other possibilities do not 
give any new information [51]. For r= 3, let us now write down the amplitude a{mi, 1712,1713) in 
the BA form given by (131). Since r =3, the sum contains 3! = 6 terms. 



a(mi,m2,m3) 



+ 
+ 
+ 
+ 



J{kimi+k2m2+k3m3)+i/2(4>i2+<l>i3+tt>23) 

„i(fcirTll+fe3TO2+fc2 m3)+i/2{4>i2+4>13+<t>32) 

^i{k2mi+kim2+k3m3)+i/2{4>2i+4>i3+4>23) 
gi{k3mi+kim2+k2m3)+i/2{4>i2+4>3i+4>32) 

„i(fe2mi+fc3m2 + feim3)+i/2(02l+03l+023) 
J(k3mi+k2m2+kim3)+i/2(<p2i+4>3i-i-4>32) 



(132) 
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For the general case n = r, i.e., r up-spins, one can proceed as in the n = 2 case and obtain the 
following equations: 

r 

e = ^(1 - coski) (133) 

i=l 

1 k- ki 

2cot-(f)j.l = COt-^ - cot— , -IT < < TT (134) 

Nki = 27rAi + ^ (135) 

j 

The Eqs. (134) are r(r-l)/2 in number. Since (/>jj = —(pij , there are only r(r-l)/2 distinct (/)'s. The 
Eqs. (135) are r in number. These r equations along with the r(r-l)/2 Eqs. of the type (134) give 
r(r+l)/2 equations in as many unknowns. Thus the equations can be expected to have solutions. 
Bethe showed that they also give the correct number of solutions. 

We now consider a particular type of solution, in which the r deviated spins form a bound 
state and move together. The following simplified analysis has been given by Ovchinnikov [52]. 

One assumes that r <<N and that only the phases 0i2,</'23, ,^r-i,r are large. One may put 

Im(j)i-i^l > without any loss of generality. Then from (135), 



Im(f>i2 = Nlmki = Nxi 
Im{(j>2i - <Ai2) = NImk2 = Nx2 

Im(pr-l,r = (136) 

One has 

Xi> for all r. Substituting in (Eq. (134)), one gets to an accuracy e~'^ 

2cot^ = 2i—-j^————: , , , = -2i (137) 

2 gi(Jie<j)i2+iim0i2) _ g— i(-Ke0i2+iim0i2) ^ ' 



Thus, 



So, 



The solution is given by 



2i = cot cot — 

2 2 

-2i = cot^ - cot^etc. (138) 



cot^ = 2i + cot%^ (139) 



coi^ = 2i/ + C (140) 
To determine C, introduce the total momentum of the deviated up-spins 

k = Y^ki = —Y^Xi 

i=l i=l 

The wave function '0 is multiplied by e**^ due to a shift — > mj + 1 and the energy levels are 
characterised by k. One gets 
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c = 

,ikl _ 

e = 



(141) 



2ri - i{e''' - 1) 

r + l{e'^- 1) 
r+ (Z- l)(e**= - 1) 

r 

y^(l - COsfc;) 
1=1 

= -(1 — cosfc) 
r 

The results can be generalised to the Hamiltonian with longitudinal anisotropy 

N 

H{J, a) = -jJ2 [S-SUi + ^(Sr Sf+i + Sf (142) 

i=l 

The anisotropy parameter a lies between and 1. The multimagnon bound states were first 
detected in the quasi-one-dimensional magnetic system CoCl2-'iH20 at pumped helium tempera- 
ture in high magnetic fields by far infrared spectroscopy [53]. Later improvements [54] made use 
of infrared HCN/DCN lasers, the high intensity of which made possible observation of even 14 
magnon bound states. 

We now describe the calculation of the ground state energy of the HAFM Hamiltonian using 
the BA. The sign of the exchange integral changes from - J to J (J > 0) and e = (^^■-^^/'^) ihe 
BA Eqns. are stih given by (133)-(135), only there is an overall 'minus' sign on the r.h.s. of (133). 
In the AFM ground state, N/2 spins are up (r=N/2) and N/2 spins down. The numbers 's are 
ordered as 

< Ai < A2 < As < (143) 

Again, Aj+i > Aj -|- 2 , for real ki 's (Case I ). There is a unique choice of the A^ 's as 

Ai = 1, A2 = 3. A3 = 5, \i^= N -1 (144) 

The ground state has total spin S = and is non-degenerate. The total ground state wave vector 
is 

N/2 

^ = E^^ = ^(l + 3 + + {N-1)) = 'Ln (145) 

If N = 4m, m positive, the wave vector is (mod 27r ); if N = 4m+2, it is tt (mod 27r). The ground 
state energy is calculated by a passage to the continuum limit: 

27 - 1 

Xj = ^ X (146) 

1 1 

1=1 







Eqs. (133)-(135) become 



2cot^<t){x, y) = cot^k{x) - cot^k{y) (148) 
1 

k(x)=2'KX+- (j){x,y)dy (149) 
2 Jo 
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N/2 

e = -Y^{l-coskj) 



N 

-— (1 — cosk{x))dx 
2 Jo 



We note that when x = y, (j) jumps from tt to - tt and 

— _ — =-2-K5{x-ij) + 2 — {cot —) 
where A =2 and B = cot^ - cot^ . Thus from (149), we get 



— = TT + cosec ( — ) — / 



dy 



4 + {cot^k{x) - cot^k{y)y 



Put 



q = cot-k 
^ 2 



and 



Po(S') = - 



dx 
dy 



Then from (152), we get the integral equation derived by Hulthen 



2 2 
-77- — = Po(5) + - 



+00 



Po{g')dg' 



, (4 + (5-ff')') 
The integral equation can be solved by the method of iterations, 

^+00 



A 



Poig')dg' 



<l+9^) J-o. (4 +(.9 -3')') 
(at the end of the calculation A is put equal to 2/n). On iteration, (156) becomes 



Po{g) = 



2 



/+00 
-00 (51 



dgi 



A 



cx, (51-5)^ + 4 

Po{g2)dg2 



.41+5?) "J-00 (4+(<?2-<7i)') 
One can continue to iterate (157) to generate higher order terms in A. 



Po{g) 



7r(l+5') 

r'+oo 



J-00 



dgi 



(4+(5i-5)^)(5? + l) 

dgidg2 



2 r 

(4 + (5i-ff)2)(4+(ff2-ffi)2)(5i + l) 
One can calculate the integrals in (158) by the method of residues. For example, 



2 /■+°° 

7^ J-00 



(4+(<?i-5)2)(5? + l) g'' + 9 
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Putting A = ^, we finally obtain the solution of the integral equation (156) as 

1 3 



IT 



+ 



(160) 



The series in (160) is uniformly convergent and one gets, 



Poig) = \sech^- (161) 



Thus, 

Po(g) 
1 + 5 

In absolute measure, the ground state energy Eg is 

NJ 



N I -r^dg = -Nln2 (162) 



Eg = — JNln2 (163) 



The low-lying excitation spectrum has been calculated by des Cloizeaux and Pearson (dCP) [55] 
by making appropriate changes in the distribution of A^s in the ground state. The spectrum is 
given by 

TT 

^ = I I' < < TT (164) 

for spin 1 states. The wave vector k is measured w.r.t. that of the ground state. A more rigorous 
calculation of the low- lying excitation spectrum has been given by Faddeev and Takhtajan [56]. 
There are S =1 as well as S = states. We give a qualitative description of the excitation spectrum, 
for details Ref. [56] should be consulted. The energy of the low-lying excited states can be written 
as E{ki, k^) = e{ki)+e{k2) with e(fcj) = ^sinki and total momentum k = ki+k2 ■ At a fixed total 
momentum k, one gets a continuum of scattering states (Figure 8). The lower boundary of the 
continuum is given by the dCP spectrum (one of the fc-s = 0). The upper boundary is obtained 
for ki = k2 = k/2 and 

k 

= TT I sin— I (165) 

The energy-momentum relations suggest that the low-lying spectrum is actually a combination of 
two elementary excitations known as spinous. The energies and the momenta of the spinous just 
add up, showing that they do not interact. A spinon is a S = 1/2 object, so on combination they 
give rise to both S=l and S=0 states. In the Heisenberg model, the spinons are only noninteracting 
in the thermodynamic limit 00. For an even number N of sites, the total spin is always an 
integer, so that the spins are always excited in pairs. The spinons can be visualised as kinks in 
the AFM order parameter (Figure 9). Due to the exchange interaction, the individual spinons get 
dclocalized into plane wave states. Inelastic neutron scattering studj^ of the linear chain S=l/2 
HAFM compound KCuF^ has confirmed the existence of unbound spinon pair excitations [57]. 
In the case of integer spin chains, the spinons are bound and the excitation spectrum consists of 
spin-w! ave-Hke modes exhibiting the Haldane gap. The BA technique described in this Section 
is the one originally proposed by Bethe. There is an algebraic version of the BA which is in wide 
use and which gives the same final results as the earlier technique. For an introduction to the 
algebraic BA method, see the Refs. [58,59]. 
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Figure Captions 

Figure 1. Ising spins on a triangular plaquette. 
Figure 2. An example of a RVB state. 
Figure 3. The Shastry- Sutherland model. 

Figure 4. Five types of interactions in the Ji — J2 — J3 — J4 — J5 model. 
Figure 5. Four columnar dimer states. 
Figure 6. A two-chain spin ladder. 

Figure 7. The continuum of scattering states (solid lines) and the bound state 
(dotted line) of two magnons. 

Figure 8. Continuum of scattering states for the S = 1/2 Heisenberg AFM chain. 
Figure 9. A two-spinon configuration in an AFM chain. 
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